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Preface 


A  preface  to  the  entire  Handbook  of  Supersonic  Aerodynamics 
appears  in  Volume  1  and  includes  a  brief  history  of  the  project.  As 
stated  in  Volume  1,  "The  primary  criterion  used  in  selecting  material 
for  the  Handbook  is  its  expected  usefulness  to  designers  of  supersonic 
vehicles.  Thus  a  collection  of  data  directly  useful  in  the  design  of 
supersonic  vehicles,  results  of  the  more  significant  experiments,  and 
outlines  of  basic  theory  are  included  ...." 

The  present  edition  of  the  Handbook,  printed  and  distributed 
by  the  Bureau  of  Ordnance,  is  being  published  in  separate  sections  as 
material  becomes  available.  The  contents  of  the  entire  work  is  given 
on  the  back  of  the  title  page.  The  majority  of  sections  of  the  Hand¬ 
book  are  now  being  prepared  by  individual  authors  for  the  Applied 
Physics  Laboratory.  Those  sections  are  noted  for  which  an  early  pub¬ 
lication  date  is  expected.  The  selection  of  material,  editing,  and 
technical  review  of  the  HandbooK  of  Supersonic  Aerodynamics  continue 
to  be  carried  out  by  an  editor  and  a  technical  reviewing  committee  at 
the  Applied  Physics  Laboratory. 

The  numbering  system  of  Section  16  is  essentially  the  same  as 
that  used  in  the  last  preceding  sections. 

Agencies  and  individuals  interested  in  the  aeronautical  sci¬ 
ences  are  invited  to  submit  and  to  recommend  material  for  inclusion 
in  the  Handbook;  full  credit  will  be  given  for  all  such  material  used. 
Regarding  the  selection  of  material  and  the  preparation  of  the  volumes 
in  the  Handbook  Series,  the  Applied  Physics  Laboratory  earnestly  soli¬ 
cits  constructive  criticisms  and  suggestions.  Correspondence  relating 
to  the  editing  of  the  Handbook  should  be  directed  to 

Editor  and  Supervisor,  Aerodynamics  Handbook  Project 

Applied  Physics  Laboratory 

The  Johns  Hopkins  University 

8621  Georgia  Avenue 

Silver  Spring,  Maryland 

Communications  concerning  distribution  of  the  Handbook  should 
be  directed  to 


Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  D.  C. 
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velocity  of  sound 
aspect  ratio 
drag  coefficient: 

lift  coefficient: 


Drae  Force 
\  P  va  S 


Lift 


Tp 


Force 
tii : 


force  coefficients  normal  to  and  parallel  to  body 
axis 

SC 

normal  force  curve-slope:  ~q£ 

specific  beat  at  constant  pressure 
characteristic  body  dimension 


E  internal  energy  flux:  — ■. x-  - 

KrJ  unit  time  x  unit  area 

Ei»  Er>  Ew  flux  of  internal  energy  on  surface  due  respectively 

to  molecules  incident  to,  and  reemitted  from  sur¬ 
face,  and  that  flux  obtainable  if  molecules  were 
reemitted  in  Maxwellian  equilibrium  with  surface 


Eint.’  Etrans. 


erf  (s  sin  8) 


portions  of  internal  energy  flux  due  respectively 
to  rotational-vibrational  energies  of  the  mole¬ 
cules,  and  that  due  to  translational  energy 


error  function: 


s  sin  6 


0 


(Values  for  this  function  are 
tabulated  in  Refs.  47  and  48.) 

erfc  (s  sin  0)  complementary  error  function: 

j^l  -  erf  (s  sin  8)j 

f  molecular  distribution  function 

h  heat  transfer  coefficient 
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,2  *2 

M'S  Mi  > 


modified  Bessel  functions  of  the  first  kind 


number  of  molecular  degrees  of  freedom 
Knudsen  number :  X/d 
thermal  conductivity 
length  along  axis  of  body 


Mach  number :  — 

a 


molecular  mass 


Ni’  \ 


number  flux  of  molecules  incident  to  and  leaving 
wall 


Nusselt  number:  -g- 

normal  stress  per  unit  surface  area;  pressure 


pi»  pr’  pw 


see  Ei?  Er,  E^  for  usage  of  subscripts 
^CP 

Prandtl  number:  — r11 


total  heat  transfer 


gas  constant 
recovery  factor : 


w(eq.) 

Tt  -  T« 


Reynolds  number : 


surface  area 

speed  ratio:  V/  V 2RT 


T  T 
w’  w(eq.) 


Stanton  number: 


S  P  V  c 


'Tw( eq. ) V 


temperature,  absolute 

wall  temperature;  adiabatic  wall  temperature  in 
the  absence  of  heat  transfer 


x-component  of  flow  velocity 

components  of  molecular  velocity  in  x<,  x«, 
directions 

velocity  of  undisturbed  flow 


xl»  x2’  x3 


local  coordinates  referred  to  a  point  on  the  sur¬ 
face  of  a  body 
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*,y,* 

a 

a 

y 

6 

0 

X 

p 

p 
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rectangular  Cartesian  coordinates 

angle  of  attack  of  body 

thermal  accommodation  coefficient: 


dE.  -  dE 
l _ r 

dE.  -  dE 
i  w 


ratio  of  specific  heats,  the  isentropic  exponent 

boundary  layer  thickness;  semi-vertex  angle  of  the 
cone 

local  angle  of  attack 
molecular  mean  free  path 
gas  viscosity 
density  (of  ambient  fluid) 


surface  reflection  coefficient  for  tangential 
momentum : 


T  -  T 
i  r 


Ti>  V  Tw 
0 
K 

General  Subscripts 

b 

t 

(eq.) 

( ideal) 

oo 


surface  reflection  coefficient  for  normal  momentum: 

pi  '  pr 
pi  ‘  pw 

tangential  stress  on  element  of  surface  area 
see  Ea,  Er,  Ew  for  usage  of  subscripts 

angle  in  cylindrical  coordinate  system 
s  sin  a 

Auxiliary  Symbols 

base 

stagnation  condition 

under  equilibrium  conditions 

theoretical  result  for  inviscid  flow 

undisturbed  free  stream  (This  subscript  is  dropped 
in  the  subsection  on  Free  Molecular  Flow.) 

condition  behind  normal  shock  wave  or  ideal  flow 
condition  at  cone  surface  behind  conical  shock 
wave 
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SECTION  16 


MECHANICS  OF  RAREFIED  GASES 


1600  Introduction 


This  section  of  the  Handbook  of  Supersonic  Aerodynamics  pre¬ 
sents  quantitative  information  of  interest  to  the  designer  working  with 
the  problems  of  flight  at  very  high  altitudes.  Together  with  Sections 
13  and  14  (Viscosity  and  Heat  Transfer  Effects)  and  Section  15  (Prop¬ 
erties  of  Gases),  the  present  section  comprises  a  single  volume.  The 
topics  proper  to  the  Mechanics  of  Rarefied  Gases  are  gathered  together 
here  into  a  single  section  of  the  Handbook  because  the  dominant  phen¬ 
omena  and  theoretical  approaches  are  different  from  those  found  in  con¬ 
tinuum  aerodynamics.  This  study  is  also  restricted  in  the  range  of 
Mach  numbers;  therefore  the  effects  of  changing  specific  heat,  dissocia¬ 
tion,  and  so  forth  (the  thermodynamics  of  which  is  treated  in  Section  15) 
are  not  considered.  Furthermore,  the  present  treatment  has  been  limited 
to  those  aspects  of  the  developing  subject  of  the  mechanics  of  rarefied 
gases  which  are  capable  of  yielding  reliable  quantitative  information. 

A  more  theoretical  consideration  of  this  subject  has  been  presented  by 
Schaaf  and  Chambre  in  an  article  entitled  "Flow  of  Rarefied  Gases”  in 
Volume  III  of  the  Princeton  Series  on  High  Speed  Aerodynamics  and  Jet 
Propulsion  (Ref.  54). 


1600.1  General  Scope  of  Contents 

It  is  possible  to  present  reliable  information  for  the  case 
of  a  gas  in  which  rarefaction  effects  are  just  becoming  evident  (slip 
flow)  and  for  the  case  of  a  completely  rarefied  gas  (free  molecular 
flow).  For  these  two  cases,  lift,  drag,  and  heat  transfer  data  for  typ¬ 
ical  bodies  are  presented  in  the  form  of  calculations  and  experimental 
results.  Also,  the  various  flow  regimes  of  rarefied  gas  dynamics  are 
discussed  and  defined. 


In  the  presentation  of  this  material  theoretical  analyses  are 
reduced  to  a  minimum  consistent  with  intelligent  use  of  the  data,  and 
reference  is  made  to  the  original  sources  which  contain  the  more  de¬ 
tailed  analyses. 
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1601.1 


1601  The  Flow  Regimes 

At  very  high  altitudes,  the  atmosphere  becomes  so  rarefied 
that  it  no  longer  behaves  like  a  continuous  fluid.  The  basic  molecu¬ 
lar  character  of  air  then  gives  rise  to  important  modifications  in 
aerodynamic  and  heat  transfer  phenomena.  Continuum  gas  dynamics  first 
must  be  modified,  then  abandoned.  The  basic  phenomena  and  theoretical 
approaches  for  highly  rarefied  flows  are  significantly  different  from 
those  for  flows  that  are  only  moderately  rarefied.  It  is  desirable 
therefore  to  divide  rarefied  gas  dynamics  into  various  flow  regimes 
(Ref.  2). 

1601.1  Qualitative  Description  of  the  Flow  Regimes 

The  characteristic  flow  in  a  highly  rarefied  gas  is  called 
"free  molecular  flow.”  In  this  regime  the  mean  free  path  is  large  com¬ 
pared  to  the  characteristic  dimension  of  an  aerodynamic  body  in  the 
flow;  and  molecules  that  impinge  on  the  body,  and  are  then  reemitted 
from  it  will,  in  general,  be  far  away  from  the  body  before  they  strike 
another  molecule.  It  follows  that  the  gas  flow  incident  on  the  body  is 
essentially  undisturbed  by  the  presence  of  the  body.  Aerodynamic  and 
heat  transfer  characteristics  depend  only  on  the  incident  flow  and  the 
average  momentum  and  energy  interaction  between  incident  molecules  and 
the  surface.  Lift,  drag,  and  heat  transfer  coefficients  may  be  calcu¬ 
lated  in  a  straightforward  way  in  terms  of  a  few  empirical  surface  in¬ 
teraction  parameters.  Experimental  results,  for  the  relatively  small 
number  of  cases  which  have  been  investigated,  are  generally  in  good 
agreement  with  these  theoretical  predictions.  However,  these  theoreti¬ 
cal  and  experimental  results  have  been  confined  to  speeds  and  tempera¬ 
tures  which  are  low  enough  to  ensure  that  no  internal  molecular  or 
atomic  energy  transformations  of  the  gas  molecules  occur  upon  striking 
the  surface.  For  very  high  velocity,  high  altitude  applications  it  is 
to  be  expected  that  molecular  dissociation,  excitation,  and  even  ioniza 
tion  will  occur  to  some  unknown  extent  on  or  perhaps  near  the  surface. 
(A  preliminary  investigation  of  free  molecular  flow  with  surface  inter¬ 
action  has  been  made  by  Sanger  (Ref.  3),  assuming  thermodynamic  equilib 
rium  at  the  surface.) 

Theoretical  methods  have  been  developed  for  correcting  free 
molecular  flow  formulae  to  account  for  the  effects  of  a  few  collisions 
(Refs.  4,  5,  and  6).  (See  Notes,  Item  1.) 

The  characteristic  flow  in  a  moderately  rarefied  gas  is 
called  "slip  flow."  This  designation  is  taken  from  the  phenomenon  of 
"slip,"  which  is  one  of  the  important  effects  known  to  occur  in  a  mod¬ 
erately  rarefied  gas  flow  and  which  is  directly  ascribable  to  its  non¬ 
continuum,  molecular  structure.  The  gas  layer  immediately  adjacent  to 
a  surface  does  not  stick  to  it,  but  instead  slips  along  it  with  a  defi¬ 
nite  velocity  proportional  to  the  product  of  the  wall  shear  stress  and 
the  molecular  mean  free  path.  A  corresponding  jump  or  discontinuity 
in  the  temperature  is  also  known  to  exist  (see  Eq.  1603.1-6). 

The  flow  regime  intermediate  between  slip  and  free  molecular 
flow  is  known  as  the  "transition  flow  regime."  It  corresponds  to  den¬ 
sities  for  which  the  mean  free  path  has  the  same  general  order  of  mag¬ 
nitude  as  the  characteristic  dimension  of  the  flow  field.  Theoretical 
treatments  are  presently  unavailable  and  even  experimental  results  (see 
Subsection  1603)  are  rather  limited  in  this  flow  regime. 
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1601.11  Real  Gas  Effects 

Flight  at  very  high  altitudes  often  involves  extremely  high 
velocities  and  resulting  high  gas  temperatures.  At  velocities  much  in 
excess  of  M  =  6  the  so-called  "real  gas"  effects  associated  with  these 
temperatures  and  due  to  the  vibration,  dissociation,  and  ionization  of 
the  gas  molecules  begin  to  be  of  importance.  The  present  section  of 
the  Handbook  is  confined  to  low  density  effects  only  and  will  not  con¬ 
sider  real  gas  effects.  Therefore  application  of  the  results  presented 
here  to  velocities  in  excess  of  M  =  6  should  be  done  with  care. 


1601.2  Significant  Parameters  for  Characterizing  the  Flow  Regimes 

Modifications  in  aerodynamic  and  heat  transfer  phenomena  due 
to  the  molecular  character  of  air  become  significant  when  the  molecular 
mean  free  path,  X,  is  comparable  to  a  relevant  characteristic  dimension, 
d,  of  the  flow  field.  The  dimensionless  parameter  called  the  Knudsen 
number,  K,  where 

K=£  ,  (1601.2-1) 


has  been  introduced  to  serve  as  a  criterion  for  determining  the  relative 
importance  of  these  rarefaction  effects.  The  Knudsen  number  is  expres¬ 
sible  in  terms  of  the  more  familiar  moduli  of  aerodynamics,  the  Mach  and 
Reynolds  numbers 


M  =  ^  and  Re  =  , 

a  ft 


(1601.2-2) 


where  V,  p,  p ,  and  a  are  the  gas  velocity,  density,  viscosity,  and  sound 
speed.  Combining  the  kinetic  theory  relation  (Ref.  1), 


M  -  |  p  A  a 


1/ST  , 

r  7r  y  ’ 


(1601.2-3) 


where  Y  is  the  isentropic  exponent,  one  obtains  a  basic  relation, 

K  =  ^  (1601.2-4) 


The  mechanics  of  rarefied  gas  dynamics  is  concerned  with  flows  in  which 
K  is  not  negligibly  small. 

The  division  of  flow  regimes  can  be  formulated  in  terms  of 
ranges  of  values  of  the  Knudsen  number  or,  by  using  Eq.  1601.2-4,  in 
ranges  of  values  of  M  and  Re. 

Free  molecular  flow  is  defined  as  that  flow  for  which  K>10. 
The  intermediate  (transition)  regime  corresponds  to  densities  for  which 
the  mean  free  path  is  of  the  same  order  of  magnitude  as  the  character¬ 
istic  dimension  of  the  flow  field. 
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Slip  Flow  is  characterized  by  a  Knudsen  number  of  a  few  per 
cent.  Here  the  relevant  characteristic  dimension  of  the  flow  field  is 
the  boundary  layer  thickness,  6,  for  large  Reynolds  numbers  or  the  body 
dimension,  d,  itself  for  small  Reynolds  numbers  (where  there  is  no 
boundary  layer) .  Since  one  has 


6 

d 


1 

Klti-  ’ 


(1601.2-5) 


at  least  two  different  Knudsen  numbers,  X/6,  or  X/d,  are  appropriate 
for  characterizing  the  slip  flow  regime.  Somewhat  arbitrarily,  but  in 
general  agreement  with  experiment,  the  slip  flow  regime  is  defined  by 
the  following  limitations: 

0.01  <  <0.1  ,  Re  >  1  ; 

/Re 

(1601.2-6) 

0.01  <  -J-  <0.1  ,  Re  <  1  . 


At  very  high  Mach  numbers  and  large  Reynolds  numbers  the  boundary  layer 
thickness  on  a  flat  plate  with  an  adiabatic  wall  behaves  more  like 


6/d  />./  — =- 
VRe 


so  it  may  possibly  be  more  appropriate  at  high  Mach  numbers  to  define 
the  slip  flow  regime  in  terms  of  the  parameter, 


1 

M  /Re 


with  perhaps  the  same  numerical  values  for  limits  as  given  above  in 
Eq.  1601.2-6.  For  the  heat  transfer  case  or  for  blunt  bodies,  still 
other  criteria  for  the  flow  regimes  are  more  pertinent  (see  Notes, 

Item  2) .  The  experimental  data  in  this  regime  are  not  yet  sufficient 
to  settle  this  question. 

These  flow  regimes,  together  with  the  equivalent  altitudes 
corresponding  to  a  characteristic  dimension  of  one  foot  are  indicated 
in  Fig.  1601.2-1.  A  detailed  presentation  of  the  variation  of  the  den¬ 
sity,  temperature,  and  composition  of  the  upper  atmosphere  with  altitude 
is  contained  in  Section  3  (General  Atmospheric  Data)  of  the  Handbook 
(see  Notes,  Item  3). 

1601.3  Surface  Interaction  Parameters 

Molecules  incident  on  the  surface  are  reemitted  from  it  with 
a  distribution  of  velocities  which  is  characteristic  of  the  surface  and 
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its  temperature,  and  the  temperature,  composition,  and  velocity  of  the 
incident  gas  flow.  For  air  incident  on  most  surfaces  of  engineering 
interest,  the  reemitted  velocity  distribution  is  almost  Maxwellian  at 
the  actual  surface  temperature.  This  makes  it  convenient  to  describe 
empirically  the  over-all  exchange  of  momentum  and  energy  in  terms  of 
departures  from  complete  Maxwellian  equilibrium.  The  "thermal  accommo¬ 
dation  coefficient,"  a,  is  defined  by 


dE,  -  dE 
i _ r 

dEi  '  dEw 


(1601.3-1) 


where  dE^  is  the  energy  flux  (per  unit  time  and  area)  incident  on  the 

surface,  dEr  the  energy  flux  reemitted  from  the  surface,  and  dE^  is 

the  energy  flux  which  would  be  reemitted  if  all  incident  molecules  were 
reemitted  in  Maxwellian  equilibrium  with  the  surface.  The  correspond¬ 
ing  coefficient, or  ,  sometimes  called  "Maxwell’s  reflection  coefficient" 
for  the  exchange  of  tangential  momentum,  is  defined  bv 


Ti  ~  Tr- 

a  -  1  r 


(1601.3-2) 


where  t\  is  the  tangential  momentum  flux  (shear  stress)  incident  on 
the  surface,  and  -r  that  for  the  reemitted  molecules.  A  third  coeffi- 

f  r 

cient,  o'  ,  to  describe  the  over-all  exchange  of  normal  momentum  is  de¬ 
fined  by 


<r’  =  Pi  ~  Pr 
^i  ^w 


(1601.3-3) 


where  p^  is  the  flux  of  normal  momentum  (pressure)  incident  on  the  sur¬ 
face,  pr  that  reemitted,  and  pw  is  the  flux  which  would  be  reemitted  if 
all  molecules  were  reemitted  in  Maxwellian  equilibrium  with  the  surface. 

The  case  in  which  5=  =  a'  =  1  is  called  "diffuse  reflec¬ 

tion,"  although  it  clearly  implies  more  than  mere  spatial  randomness  in 
the  reflection.  It  corresponds  to  an  interaction  in  which  the  incident 
molecules  achieve  complete  thermodynamic  equilibrium  with  the  surface. 
The  case  in  which  a  =  <r  =  <r*  =  0  is  called  "specular  reflection"  and 
corresponds  to  the  case  of  no  energy  or  tangential  momentum  exchange  be¬ 
tween  molecules  and  surface.  This  type  of  interaction  would  be  observed 
if  all  the  incident  molecules  were  reflected  specularly,  i.e.,  with  no 
change  in  tangential  velocity  and  simple  reversal  of  normal  velocity. 
Empirical  values  (Refs.  7  and  8)  of  a  and  <r  are  given  in  Tables  I  and 
II.  No  empirical  values  for  a  '  have  been  obtained  at  present.  It  will 
be  noted,  however,  that  for  air  incident  on  most  surfaces,  a  -  &  -  1. 

It  is  therefore  to  be  expected  that  a'  -  1  also.  (See  Notes,  Item  4.) 
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1602  Free  Molecular  Flow 

The  free  molecular  flow  regime  is  the  flow  regime  of  great 
rarefaction.  The  air  density  must  be  so  low  that  the  mean  free  path 
is  many  times  the  characteristic  dimension  of  the  body  under  considera¬ 
tion.  This  density  is  found  only  at  great  altitudes,  generally  of  the 
order  of  100  miles  or  more  (see  Fig.  1601.2-1).  The  basic  assumption 
for  theoretical  free  molecular  flow  calculations  is  that  the  flow  of 
molecules  incident  on  a  body  is  undisturbed  by  the  prese-  of  the  body, 

even  when  it  is  in  supersonic  flight.  Aerodynamic  heat  id  force 

characteristics  may  be  calculated  by  considering  separa  he  inci¬ 

dent  and  reflected  flows  of  molecules.  Experimental  agr  at  with  the 
predictions  of  theory  are  satisfactory  (Refs.  9  and  10).  For  a  brief 
discussion  of  nonuniform,  unsteady,  and  surface  interacting  free  molecu¬ 

lar  flows,  see  Schaaf  and  Chambre  (Ref.  54,  pp.  705-708). 

1602.1  Basic  Energy  Flux  Relations 

It  is  necessary  to  obtain  the  basic  energy  and  momen  im  flux 

relations  for  a  differential  element  of  area,  dS,  which  is  inclined  at 

an  angle,  0,  to  the  incident  flow.  The  coordinate  system  used  is  indi¬ 

cated  in  the  figure  below. 

x. 


*. 


X3 


It  is  assumed  that  the  incident  flow  is  in  Maxwellian  equilibrium  so 
that  its  molecular  distribution  function,  f,  is  given  by 

f  =  £  (2.BT)'3/2  exp  |  -  <u  -  V  sln  6>2  *  V  cos  9>2  +  wZ  j  , 

(1602.1-1) 

where  exp  {x}  =  e^.  Here  m  is  the  molecular  mass;  p,  T,  R,  and  V  are 

o 

respectively  density,  absolute  temperature,  gas  constant  (length  per 
2 

degree  per  time  ),  and  velocity  of  the  incident  flow;  and  u,  v,  w  are 
the  Xp  Xj,  Xg  components  of  velocity  of  a  molecule,  where  x^,  Xg, 

are  the  local  coordinates  of  the  surface.  Then  f»du»dv»dw  is  the  num¬ 
ber  of  molecules  per  unit  volume  with  velocity  components  in  the  range 
u,  u  +  du;  v,  v  +  dv;  w,  w  +  dw.  The  flux  of  energy  incident  on  the 
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surface  is  partly  of  translational  energy,  dE . ,  >,  and  partly  of 

l ( trails •  ) 

internal  energy,  dE^  ^nt  ^ .  The  former  is  given  by 


00  00 


dE 


i(  trans. ) 


j  dw  J  dv  |  |  in  (u^  +  +  w2)  f  udu  dS 


-00  -00 


=  pRT  y|T  |(s2  +  2)e"(s  sin  e>: 


+  fW  (s2  +  |)  (s  sin  0) 


1  +  erf 


(s  sin  0)j  j  dS  , 

( 1602.1-2) 


where  s  =  — — . -  is  the  molecular  speed  ratio.  The  number  flux,  dN . ,  of 

|/2RT  1 

molecules  incident  on  the  surface  is  given  by 


oo  oo  ao 

dN^  =  J"  dw  j  dv  J*  f  udu  dS 

-oc.  -oo  0 


-  ^  Y ^  |e  sin  +  yW  (s  sin  0)  ( 1  +  erf  [s  sin  0]  )  j  dS  . 

(1602.1-3) 


By  the  principle  of  equipartition  of  energy,  these  molecules  carry,  on 

the  average,  |  m  RT  units  of  internal  energy  per  molecule,  where  j  is 

the  number  of  degrees  of  freedom  for  the  internal  modes  which  partake 
in  the  energy  exchange.*  For  air  at  normal  temperatures,  j  =  2.  In 
general  j  is  related  to  the  isentropic  exponent,  /,  by 


j  = 


_  5-37 


y-  1 


The  flux  of  internal  energy  is 
_  5-3/ m  RT 


dE 


i(  int.) 


dN , 


y-  1  2 

=  o(  5  -  3 y)  (^3/2  |  e-(s  sin  0)2 

(y  -  1)  YT  *  1 

+  YW  (s  sin  0)  [1  +  erf  (s  sin  0)] 


dS 


(1602.1-4) 


(1602.1-5) 


*It  should  be  noted  that  here  each  vibrational  degree  of  freedom  must 
be  counted  as  two  ’’degrees  of  freedom"  because  potential  as  well  as 
motional  energy  can  be  stored  in  the  vibrating  system. 
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By  a  similar  line  of  reasoning,  the  total  energy  flux,  dE  ,  is  given 
by 

mRT  y  - 

dEw  =  (4  +  j)  dNw  =  gr/V-1-  T\  aRT„  dN»  .  (1602.1-6) 


w  -  s~fy'-  f)  ““‘w  W 


where  dNw  is  the  number  flux  leaving  the  surface.  For  steady  state, 


dN,  =  dN 
i  w 


(1602.1-7) 


The  net  convective  heat  transfer,  dQ,  to  the  surface,  dS,  is  given  by 


dQ  =  dE./.  v  +  dE . ,  .  .  .  -  dE  . 

i(trans.)  i(mt.)  r 


(1602.1-8) 


The  unknown  quantity,  dEr,  may  be  expressed  in  terms  of  a 
and  dEw  by  use  of  the  relation  defining  the  thermal  accommodation  co¬ 
efficient,  Eq.  1601.3-1,  to  give 


dQ  -  a  (dEi(trans-)  *  dEl(lnt))  -  5  dEw  .  (1602.1-9) 


It  will  be  noted  that  it  is  implicitly  assumed  that  a  is  the  total  ac¬ 
commodation  coefficient  for  all  modes  of  energy  which  partake  in  the 
transfer.  There  is  some  experimental  evidence  which  suggests  that 
separate  coefficients  exist  for  translational,  rotational,  and  vibra¬ 
tional  degrees  of  freedom.  The  present  treatment  will  not  consider 
such  refinements,  although  it  is  to  be  expected  that  they  will  become 
important  at  elevated  temperatures.  Expressing  the  foregoing  relation 
in  terms  of  the  appropriate  integrals, 


dQ 


_  -  — i/RT  I  ,_2  .  7  7  +  1  TW\  ( _ — ( s  sin  0)! 

-apRTy^Jis  +  -  fT"!)'  T~  V 


/¥  (s  sin  9)  £l  +  erf  (s  sin  0)J  ^  -  ■g-  e  sin 


dS 


(1602.1-10) 


This  basic  relation  may  be  integrated  over  the  surface  of  any 
convex  aerodynamic  body  to  yield  over-all  heat  transfer  characteristics. 
In  general,  the  quantities  a  and  T^  will  vary  over  the  surface  in  a  man¬ 
ner  dependent  on  the  surface  material  and  the  internal  heat  flow  charac¬ 
teristics  of  the  body.  However,  for  the  special  case  where  a  and  Tw  are 

constant  over  the  entire  surface  area,  the  indicated  integrations  have 
been  carried  out  for  a  number  of  basic  geometries  (Refs.  11  and  12). 
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The  results  which  are  of  major  interest  are  the  equilibrium  body  tem¬ 
perature,  Tw(eq  for  which  the  total  heat  transfer,  Q,  is  zero,  and 

the  heat  transfer  rate  for  body  temperatures  different  from  this. 
These  results  are  presented  most  conveniently  in  terms  of  a  thermal 
recovery  factor,  r,  and  a  Stanton  number,  St,  defined  by 


T.  -  T 
t  oo 


(1602.1-11) 


SPV  cp(Tw(eq.) 


(1602.1-12) 


where  T^  is  the  adiabatic  stagnation  temperature,  p  is  the  gas  density, 
Cp  is  the  specific  heat  at  constant  pressure,  S  is  the  total  heat  trans¬ 
fer  surface  area,  and  V  is  the  stream  velocity.  These  results  are  given 
for  the  case  of  convection  only.  It  is  to  be  noted  that  at  high  alti¬ 
tudes  the  relative  importance  of  radiative  heat  transfer  increases.  The 
complete  determination  of  thermal  characteristics  for  high  altitude 
flight  conditions  thus  must  include  a  consideration  of  radiative  heat 
transfer. 

1602.2  Basic  Momentum  Flux  Relations 

The  momentum  flux  per  unit  area  to  the  surface  element,  dS, 
is  most  conveniently  broken  up  into  a  normal  component  pi  (pressure) 

and  a  tangential  component  (shear  stress).  These  are  given  in  terms 

of  the  incident  and  reflected  molecular  flux  by 


00  00 


=  J  dw  J  dv  J  mu2  f 


-00  -00 


2  yrT  s‘ 


(s  sin  0)  e 


-(s  sin  0)‘ 


*> 

[-J  +  (s  sin  0)2J  £l  +  erf  (s  sin  0)J  > 

> 


(1602.2-1) 


00  00 


J*  dw  J"  dv  J  muvf  du 

— c,os  —  ^  e  sin  +  ylT  (s  sin  0)  £l  +  erf  (s  sin  0)J  ^  . 


2  f/Ws 


(1602.2-2) 
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1 

Using  the  reflection  coefficients,  o-  and  a  ,  one  obtains 


P  =  Pi  +  Pr 

*  (2  -  a')  p.  +  <r'  pw  (  1602.2-3) 

T  =  ri  •  Tr 

=  0-^  .  (  1602.2-4) 


The  quantity,  pw,  is  calculated  in  the  same  manner  as  the  corresponding 
energy  flux  component,  dEw: 


pw  “  7 


.  / _  dN. 

J-m  y^RT^-is1 


=  ^  >/2jtRTw  -<^e  SiD  ^  +  V ff”  (s  sin  9)  +  erf  (s  sin  0)J^ 


(1602.2-5) 


The  final  expression  for  the  normal  and  tangential  stresses  acting  on 
dS  are  thus 


p  =  f?V7F£‘8Si,ie  +  ^l/^^'<ssln  9,1 


+  (2  -  <r  )  (s^  sin^  0  +  j)  +  f  (s  sin  0)  Jjl  +  erf  (s  sin  0) 


(1602.2-6) 


T  =  — CPS-  —  ^e  s*n  +  yHF  (s  sin  0)  £l  +  erf  (s  sin  0)J  ^  , 


(1602.2-7) 


These  relations  can  be  used  (see  Refs.  13,  14,  and  15)  to  cal¬ 
culate  pressure  and  skin  friction  distributions  over  the  surface  of  any 
convex  body,  i.e.,  one  which  does  not  permit  molecular  trajectories  from 
one  surface  element  to  another.  It  will  be  observed  that  the  skin  fric¬ 
tion  depends  only  on  the  properties  of  the  incident  gas  stream  and  the 
proper  reflection  coefficient.  The  pressure,  on  the  other  hand,  depends 
also  on  the  local  body  temperature.  Thus  a  determination  of  the  thermal 
state  of  a  body  is  usually  necessary  to  determine  total  lift  and  drag 
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characteristics  in  free  molecular  flow.  Two  limiting  cases  seem  to  be 

of  most  interest.  The  first  is  the  case  of  a  "cold"  body,  taken  to  be 

at  the  ambient  g  -  temperature,  T,  and  characteristic  of  a  body  that 

is  cooled  or  is  in  the  early  part  of  its  trajectory.  The  second  is  the 

case  of  a  "hot"  body  which  has  attained  convective  thermal  equilibrium 

with  the  air  flow  and  is  at  the  temperature  T  /  *  .  Both  of  these 

w(eq.; 

cases  are  confined  to  constant  body  temperature.  However,  they  seem  to 
bracket  the  practical  possibilities  and  give  an  indication  of  the  mag¬ 
nitude  of  the  effect  of  varying  body  temperature.  For  specific  appli¬ 
cations,  it  will  probably  be  necessary  to  make  a  complete  transient 
thermal  calculation  including  both  internal  conduction  and  external 
radiation  as  well  as  convection  before  an  accurate  calculation  of  lift 
and  drag  forces  can  be  made. 

At  very  high  velocities  and  for  0> 0,  the  limiting  forms  for 
Eqs.  1602.2-6  and  1602.2-7  are  given  by 

p - —  pV2  ^(2  -  <r  )  sin2  o] 

(for  a  body  at  ambient  temperature), 


(for  a  body  at  equilibrium  temperature),  (1602.2-8) 

and 

T - -crpV2  sin  0  cos  0  .  (1602.2-9) 


These  may  be  compared  with  "Newtonian  Flow,"  an  approximation  sometimes 
used  for  a  hypersonic  continuum  flow  (Ref.  16)  in  which  it  is  assumed 
that  all  incident  normal  momentum  is  absorbed  on  the  surface,  but  that 
no  tangential  momentum  is  absorbed.  The  corresponding  expressions  are 

p( Newtonian)  *  ^  ^  6  (1602.2-10) 

T (Newtonian)  =  0  •  (1602.2-11) 

It  will  be  noted  that  the  pressure  terms  agree  for  diffuse  reflection 
for  the  "cold"  case,  otherwise  the  two  theories  are  different  and  yield 
different  results. 

For  lift  and  drag  calculation  it  is  necessary  to  resolve  the 
pressure  and  skin  friction  forces  into  the  components  normal  and  paral¬ 
lel  to  the  direction  of  V.  When  the  differential  lift  and  drag  forces 
acting  on  dS  are  denoted  by  dFL  and  dFp  respectively,  the  proper  rela¬ 
tions  are 


dFL  =  (p  cos  0  -  r  sin  0)  dS  (1602.2-12) 

dFjj  =  (  t  cos  0  +  p  sin  0)  dS  . 


(1602.2-13) 
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The  total  lift  and  drag  forces  are  then  obtained  by  integration  over 
the  surface  area.  It  is  convenient  to  express  these  results  in  terms 
of  lift  and  drag  coefficients,  and  CD,  defined  by 


C,  = 


C„  = 


Lift  Force 

S  pv2s 


Drag  Force 

*pv»s 


(1602.2-14) 


(1602.2-15) 


where  S  is  the  reference  area.  In  keeping  with  usual  aerodynamic  ter¬ 
minology  the  area,  S,  has  been  taken  as  the  maximum  cross  sectional 
area  of  the  aerodynamic  configuration  normal  to  the  body  axis,  except 
as  otherwise  noted. 


1602.3  Heat  Transfer  for  Typical  Bodies  ( Calculated) 

The  results  of  integrating  Eq.  1602.1-10  over  the  total  sur¬ 
face  for  several  different  geometries  are  presented  below. 

1602.31  The  Flat  Plate  at  Angle  of  Attack,  a 

The  results  for  the  flat  plate  with  front  and  back  surfaces 
insulated  from  one  another  are,  for  the  front  surface  only. 


r  = 


(y+  l)si 


2s2  +  1 


2 

1  +  Vn  (s  sin  a)  [  1  +  erf  (s  sin  a)]  e^s  sin  c ^ 


(1602.31-1) 


St 


(e"<s  sln  “>2  +  ^(s  sin  a) 
4  \fv  sy  ^ 


[1  +  erf  (s  sin 


(1602.31-2) 

These  results  are  plotted  in  Figs.  1602.31-1 
and  1602.31-4.  The  formulas  for  the  back 
side  are  obtained  by  replacing  a  and  -a  and 
are  plotted  in  Figs.  1602.31-2  and  1602.31-5. 
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The  results  for  the  flat  plate  with  front  and  back  surfaces 
in  thermal  contact  are: 


r  = 


( r  *  Ds‘ 


2s*  +  1  - 


A. 


1  +  YtT  (s  sin  a)  erf  (s  sin  a)  e^s  sin  J 

( 1602.31-3) 


St 


=  ^  i  e  ^ s  sin  +  y'W  s  sin  Of  erf  (s  sin  a)  \  . 

4  YtT  sY  L  J 


( 1602.31-4) 


These  results  are  plotted  in  Figs.  1602.31-3 
and  1602.31-6. 


1602.32  The  Infinite  Right  Circular  Cylinder  at  Angle  of  Attack,  a 


r  = 


(Y  +  l)s‘ 


[f2  +  2s2(1  *f2)]  I0(^)  +  [{2(1  +  2s2)]  I,(^) 

(1  +  i2)  Io<T*  +  C2  *1 


St 


4  y  yT  s 


(l 


*«2> 


<^r> 


( 1602.32-1) 


(1602.32-2) 


p2  f2 

where  £=  s  sin  a  and  Iq(  2^  and 

are  modified  Bessel  Functions  of  the  first 

kind,  i.e.,  IQ(x)  =  i  nJfl( ix)  =  inJn(-ix). 

Values  for  these  functions  of  order  0  and 
1  are  tabulated  in  Ref.  49.  The  total 
heat  transfer  area  is  taken  as  the  curved 
surface  area,  TrLd,  where  L  is  the  length 
of  the  cylinder. 


For  a  -  0,  these  results  reduce  to: 


r  ( a  =  0) 


_  2  Y 

Y  +  1 


=  (r  +  D« 

4  y  ylT  s 


St  (a  =  0) 
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For  purposes  of  graphical  representation  it  is  convenient  to  let 


and 


where 


) 


and 


^-y— -  r  =  2  +  sin^  a  *  h  ((  ) 


7  _  sin  ot  . 

*  1}  ~  4  y¥ 


st  =  •  g  (O  , 


,2  c  2 

,  _  l0  (V  +  X1  (V 

h(4) - 75 - 7T 

(1  +  e2)  lfl($  +  1^) 


=  [(|-  +  £)  Vt*  +  £  ] 


e" 


( 1602.32-3) 


( 1602.32-4) 


( 1602.32-5) 


(1602.32-6) 


See  Figs.  1602.32-1,  1602.32-2,  and  1602.32-3. 


1602.33  The  Sphere 


r 


r  = 


y  v 


(2s^  +  1)  (  1  +  g  ierfc(s) j  +  — — g— I  erf(s) 


^1  +  i  ierfc(s)j  +  erf(s) 


St  = 


8 


^  {s*  +  s  ierfc(s)  +  ^  erf(s)  ^  , 


(1602.33-1) 

(1602.33-2) 


where  ierfc(s)  is  the  complementary 
error  function  integrated  over  the 

2 

total  heat  transfer  area,  wd  .  See 
Figs.  1602.32-1  and  1602.32-2. 


NAVORD  REPORT  1488  (Volume  5) 

1602.34  Mechanics  of  Rarefied  Gases  February  1959 


1602.34  The  Cone  at  Angle  of  Attack,  a  ( Conical  Surface  Only) 


where  the  height  of  the  cone  is  L  and  the  semi-vertex  angle  is  6.  The 
total  heat  transfer  area  is  the  curved  surface  area,  and 


sin  0  -  -  cos  6  cos  0  sin  a  +  cos  a  sin  6  . 


(1602.34-2) 


The  Stanton  number,  St,  for  the  cone  may  be  obtained  by  substituting 
Eq.  1602.34-1  into  Eq.  1602.1-12.  The  value  of  r  will  be  found  as  be¬ 
fore  by  equating  the  integral  in  Eq.  1602.34-1  above  to  zero,  i.e., 

•p 

Q  =  0,  and  by  making  the  following  substitution  for  : 


T_w(e<*«) 

T 


1  +  r 


LI  -  11  2 


(The  cone  at  zero  angle  of  attack  may  be  treated  by  using  the  flat  plate 
result  with  a  replaced  by  6,  see  Eqs„  1602.31-1  through  1602.31-4.) 


1602.35  Practical  Application  of  the  Heat  Transfer  Equations 

Composite  bodies  that  are  convex  may  be  treated  in  relatively 

simple  fashion.  If  the  thermal  characteristics,  r.,  and  St.,  and  areas 

J  3 

Sj,  where  j  =  1,  ...n,  of  a  series  of  n  basic  shapes  are  known,  the 

heat  transfer  characteristics  for  the  combined  configuration  can  be  ob¬ 
tained  from 


n 

r  =  sk  £  r j  sj  st j 

*3  • 

n 

*  S  £  St  s  , 
j  =  1  J  J 


St 
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provided  the  temperature  of  the  composite  body  is  constant.  The  areas, 


posite  shape  exposed  to  the  flow. 


As  indicated  in  Subsection  1602.1  Page  4,  the  values  of  r  and 
St  in  the  equations  above  are  average  or  over-all  values  of  these  param¬ 
eters.  In  order  to  calculate  the  heat  transfer  or  temperature  at  speci¬ 
fic  points  on  the  configuration,  it  will  be  necessary  to  integrate 
Eq.  1602.1-10  for  specific  values  of  0. 

Also,  it  will  be  observed  that  the  thermal  recovery  factors 
are  in  general  greater  than  unity.  Body  equilibrium  temperatures,  with 
respect  to  convection,  are  greater  than  the  adiabatic  stagnation  tem¬ 
perature  of  the  gas  flow.  This  striking  phenomena  is  one  of  the  more 
important  and  unusual  features  of  the  free  molecular  flow  regime.  It 
will  also  be  observed  that  heat  transfer  to  a  ’’cold"  body  is  propor- 

3 

tional  to  pV  S  for  high  velocity  flight.  This  is  in  accord  with  meteor 
observations. 


1602.4  Aerodynamic  Forces  for  Typical  Bodies  ( Calculated) 

The  results  of  integrating  Eqs.  1602.2-12  and  1602.2-13  over 
the  total  surface  area  of  several  different  typical  geometrical  shapes 
are  presented  below.  These  calculations  have  been  carried  out  under 
the  assumption  that  T  ,  <r  ,  and  <r'  are  constant  over  the  entire  surface. 


1602.41  The  Flat  Plate  at  Angle  of  Attack,  a,  with  Both  Sides  Ex¬ 
posed  to  th^  Flow 

C,  =  2(2  r  a-  -  cos  „  sln  a  e~(s  sin  or)2  +  £  l/E!>  „  cos  „ 

L  yTTs  s  T 


+  COS 


a  2(2  - 


*  V  .2  2  -  O' 

cr  ~  <r  )  sin  a  +  - x — 


erf(s  sin  a) 

(1602.41-1) 


CD  = 


[~2(  2  -  <r* 

L  yir  s 


sin^  a  +  ■  -  —  cos' 
yn  s 


2  a  j  e'(s  si”  “)2  +  -f-  V^sin2  a 


a  ( 


t  o  1  2 

+  2  sin  a  ( 2  -  or  )  (sin  a  +  — y)  +  a  cos  a 


2s 


erf(s  sin  a) 


(1602.41-2) 


The  reference  area  is  taken  as  the  area 
of  one  side  of  the  plate.  These  results 
are  plotted  in  Figs.  1602.4-1  through 
1602.4-6. 
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1602.42  The  Sphere 


2-0- 

2? 


+  <r  f  4s^  +  4s^  -  1 

-  | - 


er,(s,  ♦  l?-s2  *  1)  e-'2 


7  s  f  T 


(1602.42-1) 


The  reference  area  is  the  projected  area. 
This  equation  is  plotted  in  Fig.  1602.4-7. 


1602.43  The  Cylinder  at  Angle  of  Attack,  a 


i  3/0  2 

■Z_i - «a_s  ♦  2i£_s  ( „  +  4  -  2  ,• ) 


r  TV 


+ 1  t2)  i0(^)  +  +  i  e2)  } 


(1602.43-1) 


£ 

;^T_e  Y  cos  a  |(1  +  e2}  1^)  +  (2  l^)  j  >  ( 1602.43-2) 


where  the  reference  area  is  diameter,  d, 
n  times  unit  length. 

-—cD 


Here, 


CN  =  CL  COS  a  +  CD  Sin  a  ’  Ca  =  CD  COS  a  ~  CL  Sin  a 


(1602.43-3) 


and  the  notation  is  otherwise  the  same  as  for  the  cylinder  heat  trans¬ 
fer  case.  For  purposes  of  graphical  representation,  it  is  convenient 
to  let 


cA  =  VT  tr  cos  a  sin  a  g(  £ ) 


(1602.43-4) 
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;n  -  -If-  "%n  ^T-  *  (4  +  <r  -  2  <r’)  sin2  .  !((|  , 


(1602.43-5) 


where 

f(c )  =  v*~  +  i  c  >  +  +  4>  V' ]  e  * 72 

(1602.43-6) 

and  g(£)  is  defined  in  Eq.  1602.32-6.  See  Figs.  1602.32-3  and  1602.4-7, 
1602.44  The  Cone  at  Angle  of  Attack,  a 


Cv  =  — « r -  (p  cos(n,x)  +  rcos(t,x))  d0 

N  pV^ff  sin  6 


(1602.44-1) 


CA  = 


pv  tt  sin  5 


ii 

|  (p  cos(n,y)  +  rcos( t, y) )  d0  (1602.44-2) 


where  6  is  semi-vertex  angle,  p  and  r  are  the  expressions  given  in 
Eqs.  1602.2-6  and  1602.2-7,  and  sin  &  in  these  expressions  is  the  same 
as  that  given  in  Eq.  1602.34-2.  The  coefficients  CN  and  are  referred 

to  the  base  area.  The  angles  and  coordinate  systems  used  in  the  cone 
equations  are  indicated  in  the  figure  below. 
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The  direction  cosines  are: 


cos(n,x)  =  cos  6  cos  0 


( 1602.44-3) 


cos(n,y)  =  sin  6 


(1602.44-4) 


cos 


2  2 

I  .  *  sin6  (sing  sinfi  +  cosa  cosfi  cos0)  +  sing  cos  6  sin  0 

^t,x;  .  .  —  —  ...  — .  —  —  —  ...  — 


i  Vc 


.2*  _,_2 


cos  6  sin  0  +  (sin6  sina  +  cos6  cos0  cosa)' 


(1602.44-5) 


cos(t,y)  = 


2  2 

cosa  cos  6  sin  0  +  cos6  cos0  (cosa  cos6  cos0  +  sina  sinC) 

— 

+  y  cos  6  sin  0  +  (sinfi  sina  +  cos6  cos0  cosa) 


(1602.44-6) 


The  values  of  cos(n,y)  and  cos(t,y)  are  always  positive.  The  value  of 
cos( n, x)  goes  from  positive  to  negative  values  as  0  goes  from  0  to  v 
while  cos(t,x)  goes  from  negative  to  positive  values  as  0  goes  from 
0  to  n.  The  +  sign  in  Eqs.  1602.44-5  and  1602.44-6  is  chosen  so  that 
this  is  the  case.  These  results  can  be  expressed  in  terms  of  CL  and  Cp 

using  Eq.  1602.43-3  and  are  so  plotted  in  Figs.  1602.4-8  through 
1602.4-11.  (See  Notes,  Item  5.) 


1602.45  Example  of  £  Composite  Body 

Composite  bodies  which  are  convex  may  be  treated  similarly  to 
the  heat  transfer  case.  If  a  configuration  consists  of  n  different 
parts,  each  of  which  has  a  known  lift  and  drag  coefficient, 

(Cjj)j,  j  =  1,  ...n,  then  the  lift  and  drag  coefficients,  CL  and  CD,  of 

the  total  are  given  by 

n 

CL  =  5  J?,  (cl)j  sj  (1602.45-1) 


n 


where  S.  and  S  are  all  the  appropriate  reference  areas. 


(1602.45-2) 
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1603  Slip  Flow 

The  slip  flow  regime  is  the  regime  of  moderate  rarefaction, 
corresponding  to  densities  at  which  the  gas  first  begins  to  exhibit  its 
non-continuum  molecular  structure.  The  mean  free  path  is  a  few  per  cent 
of  the  characteristic  dimension  of  the  flow  field.  Typical  altitudes  of 
interest  are  in  the  range  from  20-50  miles,  depending  on  the  Mach  number 
and  size  under  consideration. 


1603.1  The  Velocity  Slip  and  Temperature  Jump  at  Body  Surfaces 

The  slip  flow  regime  takes  its  name  from  the  phenomenon  of 
"slip"  according  to  which  a  gas  at  low  pressures  does  not  stick  to  a 
surface  but  rather  slips  along  it  at  a  definite  rate  determined  by  the 
mean  free  path,  the  shear  stress,  and  other  quantities. 

A  very  simple  analysis  will  serve  to  indicate  how  this  comes 
about.  Consider  a  gas  moving  with  velocity  U(y)  in  the  x-direction 
past  a  surface  lying  on  the  x-axis.  The  layer  of  gas  immediately  adja¬ 
cent  to  the  surface  consists  of  molecules  of  which  approximately  half 
have  just  been  emitted  from  the  surface,  while  the  remainder  have,  on 
the  average,  come  from  a  layer  a  mean  free  path  distant  from  the  sur¬ 
face.  The  average  tangential  velocity  of  these  incident  molecules  will 
be  given  by 


U(  A.)  =  U(0).  +  A 


dU 

dy 


(1603.1-1) 


If  cr  is  the  reflection  coefficient  for  tangential  momentum,  the  average 
tangential  velocity  component  of  the  reemitted  molecules  is  given  by 

(1  -  <r  )U(A)  =  (1  -  <r  )  [u(0)  +  A  (1603.1-2) 

The  velocity  of  the  gas  layer  itself  is  thus 


U(0)  =  £  }  U(0)  +  A  +  (1  -  <r  ) 


From  this  it  follows  that 

U(0)  =  A  ^  .  (1603.1-4) 


given  by 

[««♦>*]}• 

(1603.1-3) 


At  normal  pressures,  A  is  so  small  that  U( 0)  =  0  to  a  very  good  approxi¬ 
mation.  At  very  low  pressures,  however,  A  becomes  large  enough  to  make 
possible  an  appreciable  slipping  velocity.  It  will  be  noted  that  the 
phenomenon  occurs  even  for  completely  diffuse  reflection,  i.e,,  <^=1. 

By  similar  analyses,  the  velocity  slip  due  to  a  tangential 
temperature  gradient  and  the  discontinuity  between  gas  and  surface 


532927  0-59-3 


NAVORD  REPORT  1488  (Volume  5) 
1603.1  Page  2  Mechani 


Mechanics  of  Rarefied  Gases 


February  1959 


temperature  for  the  case  of  heat  flow  may  be  obtained.  The  complete 
expressions  (Ref.  24)  are 


U(gas)  =  ^2-  X  (§£)  +  — Lr 

a  ay  w  Vzrf 


[«  si  <ln  T>] 


T( gas)  -  T( wall)  = 


2-5  .  2  X  _X_  (Ql\ 

&  y  +  l  pr  'dy  , 


( 1603.1-5) 


(1603.1-6) 


where  a  is  the  thermal  accommodation  coefficient,  Pr  is  the  Prandtl  num¬ 
ber,  Pr  =  and  y  is  the  isentropic  exponent.  Special  attention 

should  be  called  to  the  second  term  on  the  right  hand  side  of 
Eq.  1603.1-5,  according  to  which  a  flow  is  induced  along  a  surface  in 
the  direction  of  increasing  temperature.  This  is  called  "thermal  creep" 
and  can  be  of  considerable  importance  in  very  low  pressure  instrumenta¬ 
tion.  These  relations  are  approximate  in  the  sense  that  they  apply  only 
for  moderate  departures  from  equilibrium.  If  either  the  velocity  or 
temperature  gradients  are  very  large,  more  complicated  relations  appar¬ 
ently  exist.  In  the  slip  flow  range,  as  defined  in  Subsection  1601.2, 
however,  they  are  quite  satisfactory  and  have  been  well  verified  by  nu¬ 
merous  experiments.  It  will  be  observed  that  slip  and  temperature  jump 
phenomena — by  themselves — have  the  effect  of  reducing  skin  friction  and 
heat  transfer,  acting  in  many  ways  as  if  an  extra  film  of  air  of  the 
order  of  a  mean  free  path  thick  rested  on  the  solid  surface. 


1603.2 


Differential  Equations  for  Slip  Flow 


There  is  still  some  question  as  to  the  appropriate  set  of  dif¬ 
ferential  equations  to  describe  a  gas  flow  in  the  slip  flow  regime.  The 
most  recent  experimental  and  theoretical  evidence  seems  to  indicate  that 
the  Navier-Stokes  equations,  together #with  a  bulk  viscosity  modulus  of 
the  order  of  two-thirds  of  the  shear  modulus,  are  satisfactory  for  air 
in  the  slip  flow  regime.  Other,  more  complicated  sets  of  equations, 
known  as  the  Burnett  equations  (Refs.  17  and  18),  the  Thirteen  Moment 
equations  (Ref.  19),  and  others  (Ref.  20)  have  been  formulated  for 
monatomic  gases  by  means  of  the  kinetic  theory  of  gases.  However,  these 
sets  of  equations  do  not  appear  to  be  as  yet  adequately  substantiated  by 
either  experimental  evidence  (Refs.  21  and  22)  or  by  theoretical  evi¬ 
dence  (Ref.  23)  for  monatomic  gases,  let  alone  air.  They  are  in  general 
not  superior  to  the  simpler  Navier-Stokes  relations.  (See  Notes,  Item  6. 
For  a  more  detailed  discussion  of  the  Thirteen  Moment  and  the  Burnett 
equations  and  some  of  the  difficulties  associated  with  these  equations, 
see  Schaaf  and  Chambre,  Ref.  54,  pp.  710-718.)  In  the  transition  regime 
no  set  of  macroscopic  differential  equations  yet  proposed  seems  to  be 
valid. 


1603.3  Viscous  Interact  ion  Effects — Hypersonic  Flow  Results 

From  the  basic  relationships  indicated  in  Eq.  1601.2-6,  i.e., 

X  M 

the  limitation  on  K  =  -  , it  follows  that  non-continuum  phenomena 

0  Vte 
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will  occur  in  general  at  high  Mach  numbers  or  at  low  Reynolds  numbers, 
or  perhaps  both.  But,  to  some  extent,  this  is  also  the  same  range  of 
parameters  for  which  "interaction  effects"  between  the  viscous  and  in- 
viscid  parts  of  a  flow  field  are  known  to  occur.  Generally  speaking, 
boundary  layers  are  thick;  they  force  the  inviscid  flow  past  the  body 
outwards  a  distance  the  order  of  the  displacement  thickness,  thus  pre¬ 
senting  an  effectively  thicker  body  to  the  flow;  and  this  perturbs  the 
surface  pressure  distribution  which  in  turn  feeds  back  into  the  boundary 
layer  flow  to  increase  the  skin  friction.  These  are  really  continuum 
flow  effects.  The  effects  of  slip  and  temperature  jump  at  the  surface 
are  non-continuum  phenomena  but  are  of  importance  for  aerodynamic  appli¬ 
cations  usually  only  in  the  presence  of  strong  interaction  effects. 

Slip  and  temperature  jump  have  the  effect  of  reducing  skin  friction 
and  heat  transfer,  while  interaction  phenomena  have  the  effect  of  in¬ 
creasing  skin  friction  and  heat  transfer.  Thus,  there  is  a  competition 
of  two  tendencies  of  opposite  sign.  No  complete  theoretical  solution 
for  a  flow  situation  of  aerodynamic  interest  is  presently  available. 
However,  a  growing  body  of  empirical  data  relating  to  heat  transfer  and 
aerodynamic  characteristics  for  basic  geometrical  shapes  in  slip  and 
transition  flow  does  exist  and  seems  to  indicate  that,  with  decreasing 
density,  interaction  effects  are  initially  dominant  but,  at  lower  den¬ 
sities,  the  trend  is  reversed  and  slip  and  temperature  jump  effects 
dominate.  This  body  of  data  will  be  summarized  below. 


1603.4  Heat  Transfer  for  Typical  Bodies  ( Experimental) 

Heat  transfer  characteristics  are  most  conveniently  presented 
in  terms  of  a  thermal  recovery  factor,  r,  and  a  Nusselt  number,  Nu: 

r  =  (1603.4-1) 

At  00 

Nu  =  ■*—  ,  (1603.4-2) 

where  h  is  given  from 

«  *  hS  <T.(e,.)  -  T„>  ;  (1603.4-3) 

and  where  k  is  the  thermal  conductivity  of  the  gas,  T,^  is  the  free 

stream  gas  temperature,  T^  is  the  stagnation  temperature,  T^  is  the 

wall  temperature,  T^ ^  is  the  wall  temperature  under  equilibrium 

conditions,  h  is  the  heat  transfer  coefficient,  and  d  is  a  characteris¬ 
tic  dimension.  The  body  geometries  for  which  values  of  these  param¬ 
eters  have  been  determined  are  denoted  in  the  following  subsections. 


1603.41  The  Sphere 

Heat  transfer  measurements  have  been  obtained  for  spheres 
(Refs.  25  and  26)  over  the  ranges  2.7<M<6,  15<Re<7000  and  0<M<0.7, 
2<Re<100.  These  results  are  presented  in  Figs.  1603.41-1,  1603.41-2, 
and  1603.41-3. 
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In  the  reduction  of  these  data,  free  stream  properties  are 
used  for  the  evaluation  of  the  Nusselt  and  Reynolds  numbers  for  the 
subsonic  data  and  for  the  Reynolds  number  for  the  supersonic  recovery 
factor  data;  properties  and  flow  conditions  behind  a  normal  shock  wave 
are  used  for  the  evaluation  of  the  Nusselt  and  Reynolds  numbers  for  the 
supersonic  heat  transfer  data;  and  the  total  surface  area  of  the  sphere, 
2 

JTd  ,  is  taken  as  the  reference  area. 


1603.42  The  Right  Circular  Cylinder 

Heat  transfer  results  have  been  obtained  for  right  circular 
cylinders  cross  stream  to  the  flow  (Refs.  27,  28,  and  29)  over  the 
range  1.1<M<4.5,  0.1<Re2<  200.  These  results  are  presented  in 

Figs.  1603.42-1,  1603.42-2,  and  1603.42-3. 

The  Reynolds  number,  Reg,  is  based  on  the  velocity,  density, 

and  viscosity  behind  a  normal  shock  wave,  and  the  characteristic  dimen¬ 
sion  is  the  cylinder  diameter.  The  reference  area  (per  unit  length  of 
the  cylinder)  is  taken  as  ;rd.  (See  Notes,  Item  7.) 


1603.43  The  Cone 

Heat  transfer  results  have  been  obtained  for  cones  at  zero 
angle  of  attack  (Ref.  30)  over  range  10°<6<60°,  2.16<M<3.54, 
80<Re<1800  (where  6  is  the  semi-vertex  angle).  These  results  are 
presented  in  Figs.  1603.43-1  and  1603.43-2. 

The  Reynolds  number  is  based  on  the  free  stream  velocity, 
density,  and  viscosity,  and  the  characteristic  dimension  is  the  cone 
slant  length.  The  reference  area  is  the  total  conical  surface  area. 

The  results  for  cones  were  obtained  from  experiments  using 
an  insulated  base  area  shielded  from  the  flow;  therefore,  they  apply 
to  the  conical  surface  portion  only. 


1603.5  Aerodynamic  Forces  for  Typical  Bodies  ( Experimental) 

The  lift  and  drag  characteristics,  as  well  as  some  surface 
pressure  distributions,  which  have  been  determined  for  a  number  of 
basic  geometrical  shapes  in  the  slip  flow  and  transition  flow  regimes 
are  denoted  and  briefly  discussed  in  the  following  subsections. 


1603.51  The  Sphere 

Drag  characteristics  have  been  determined  for  spheres 
(Refs.  31  and  32)  over  the  ranges  M^O,  20<Re<300  and  M~2.5, 
8<Re2<  400.  These  results  are  presented  in  Fig.  1603.51-1  in  terms 

of  a  drag  coefficient,  Cn,  where  the  reference  area  is  the  projected 
2  u 

area,  — .  The  Reynolds  number,  Re,  is  based  on  the  velocity,  den¬ 
sity,  and  viscosity  behind  a  normal  shock  wave  and  the  diameter,  d, 
of  the  sphere. 
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1603.52  The  Flat  Plate 

Drag  forces  on  flat  plates  at  zero  angle  of  attack  (Refs.  33 
and  34)  have  been  determined  over  the  range  0<M<3.80,  3  <  Re  <  2000. 
These  results  are  presented  in  Fig.  1603.52-1  in  terms  of  a  drag  coef¬ 
ficient  based  on  total  drag  force,  free  stream  flow  conditions,  and  a 
reference  area  taken  to  be  one  side  of  the  flat  plate.  The  Reynolds 
number  is  based  on  the  free  stream  flow  and  the  length  of  the  plate. 
These  results  apply  to  flat  plates  of  finite  length,  and  thus  contain 
the  induced  drag  due  to  the  wake  at  the  trailing  edge  (Ref.  35).  Span 
and  thickness  effects,  however,  have  been  essentially  eliminated,  so 
that  the  results  correspond  to  infinitely  thin  plates  with  infinite 
aspect  ratio. 

Lift  forces  for  flat  plates  (Ref.  36)  at  an  angle  of  attack, 
at,  have  been  obtained  over  the  range  3.70  <M<  4.13,  630 <  Re <2400. 
These  results  are  given  in  terms  of  a  normal  force  coefficient,  C^, 
given  by 


„  _  Normal  Force 

CN - I  2 - 

^pV  S 


(  1603.52-1) 


where  S  is  the  plate  area,  and  p  and  V  are  the  free  stream  density  and 
velocity  respectively.  The  Mach  number  and  Reynolds  number  are  based 
on  free  stream  flow  conditions,  and  the  characteristic  dimension  is  the 
length  of  the  flat  plate.  It  was  observed  experimentally  that  C„  varied 
linearly  with  a  in  the  range  0  <  a  <  8°,  so  that  the  quantity 


C  =  ~-  (1603.52-2) 
a 

could  be  determined.  For  the  inviscid  supersonic  case  and  for  flat 
plates  of  finite  aspect  ratio,  A,  this  quantity  can  be  calculated  on 
the  basis  of  linearized  theory,  i.e., 


CN  ( Lin.Th. ) 


The  experimental  results  are  hence  presented  in  Fig.  1603.52-2  in  terms 
of  the  ratio  /C^  (Lin.Th.).  (See  Notes,  Item  8.) 


1603.53  The  Cone 

Cone  surface  pressure  variation,  due  to  the  thick  laminar 
boundary  layer,  has  been  determined  (Refs.  37  and  38)  over  the  range 
3.7<M<5.8,  185  <  Re  <  30,  000.  The  cones  had  3°  and  5°  semi-vertex 
angles  and  were  at  zero  angle  of  attack.  The  local  Reynolds  number 
was  based  on  the  flow  behind  an  ideal  conical  shock  and  the  cone  slant 
length  from  vertex  to  pressure  orifice.  These  results  are  presented 
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in  Fig.  1603.53-1  in  terms  of  the  variation  of  a  local  fractional  pres¬ 
sure  increase  with  the  local  Reynolds  number.  The  fractional  pressure 
increase  is  expressed  as 


( ideal) 


P( ideal) 


( 1603.53-1) 


where  p  is  the  cone  surface  pressure  and  p 
inviscid  pressure  on  the  cone  surface. 


( ideal) 


is  the  theoretical 


Cone  drag  forces  for  cones  at  zero  angle  of  attack  and  15° 
semi-vertex  angle  have  been  determined  ( R«  i  .  39)  over  the  range  M~2, 
75<Re<  1800;  M~  4,  500<  Re  <7000.  Only  the  tip  region  was  investi¬ 
gated;  hence  the  shoulder  effect  was  not  present  and  the  cones  were 
essentially  infinite  in  length.  The  results  are  presented  in 
Fig.  1603.53-2  in  terms  of  the  increase  in  the  drag  coefficient, 

Cp  -  ideai)»  over  and  above  the  inviscid  supersonic  drag 

c  lent,  ^  i  \>  at  ttie  test  Mach  number.  The  drag  coefficient 


coef f i- 

C„, 


is 


D(  ideal)’ 

based  on  free  stream  flow  conditions,  the  projected  cone  surface  area, 
and  the  total  measured  drag.  It  thus  consists  partly  in  skin  friction 
drag  and  partly  in  wave  drag.  The  curves  indicated  in  Fig.  1603.53-2 
are  taken  from  theory  (Ref.  40).  (See  Notes,  Item  9.) 


1603.54  Base  Pressures 

The  variation  of  base  pressure  on  a  cone-cylinder  configura¬ 
tion  at  zero  angle  of  attack  has  been  determined  (Refs.  41,  42,  and  43) 
over  a  very  large  range  of  Reynolds  numbers.  The  Reynolds  numbers  were 
based  on  free  stream  flow  conditions  and  the  model  length.  The  con¬ 
figuration  consisted  of  a  cone  of  30°  semi-vertex  angle  followed  by  a 
cylindrical  afterbody  so  that  the  over-all  length  diameter  ratio  was 
10:3.  The  results  are  presented  in  Fig.  1603.54-1  in  terms  of  the 
ratio  of  base  pressure,  pb,  to  ambient  pressure,  p^,  versus  the  Reynolds 

number.  In  the  lower  Reynolds  number  range  the  pressure  varied  con¬ 
siderably  with  radius  over  the  base  region.  The  base  pressures  in  the 
figure  are  area  averages. 


1603.55  Impact  Probe  Measurements 

At  low  Reynolds  numbers  the  pressure  at  the  forward  stagna¬ 
tion  point  of  a  two-  or  three-dimensional  body  increases  above  the  value 
of  the  local  stagnation  pressure.  The  effect  depends  upon  the  Mach  num¬ 
ber  and  Reynolds  number,  as  well  as  the  geometry  of  the  body.  This  ef¬ 
fect  is  particularly  important  in  connection  with  the  interpretation  of 
impact  or  total  head  probes.  Early  experiments  in  incompressible  flow 
were  carried  out  by  Homann  (Ref.  44).  More  recent  experiments  have 
been  carried  out  by  Kane  and  Maslach  (Ref.  45)  and  Sherman  (Ref.  46). 

The  results  (Ref.  46)  for  subsonic  flow  are  presented  in  Fig.  1603.55-1 
in  terms  of  the  pressure  coefficient  versus  Reynolds  number,  based  on 
free  stream  flow  conditions  and  the  probe  diameter.  The  results  for 
the  supersonic  case  (Ref.  46)  are  presented  in  Figs.  1603.55-2  and 
1603.55-3  in  terms  of  the  ratio  of  the  measured  impact  pressure  to  the 
ideal  value,  i.e.,  the  stagnation  pressure  behind  a  normal  shock  wave 
at  the  free  stream  Mach  number.  (See  Notes,  Item  10.) 
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SECTION  16  MECHANICS  OF  RAREFIED  GASES 


Notes 


Some  of  the  more  important  theoretical  and  experimental  re¬ 
sults  which  have  appeared  since  the  preparation  of  this  manuscript  and 
which  could  not  be  included  in  the  body  of  the  text  are  listed  below. 

1.  Calculations  of  the  drag  and  heat  transfer  for  a  sphere  in  the  near- 
free  molecular  flow  range  are  given  in 

Baker,  R.  M.  L.  and  Charwat,  A.  F.  "Transitional  Cor¬ 
rection  to  the  Drag  of  a  Sphere  in  Free  Molecule  Flow," 

Phys.  of  Fluids.  Vol.  1  (March-April  1958),  pp.  73-81. 

Hammerling,  P.  and  Kivel,  B.  "Heat  Transfer  to  a 
Sphere  at  the  Transition  from  Free  Molecule  Flow," 

Phys.  of  Fluids.  Vol.  1  (July-August  1958),  p.  357. 

2.  A  discussion  of  the  regimes  of  flow  for  highly  cooled  blunt  bodies 
is  given  in 

Adams,  Mac  C.  and  Probstein,  R.  F.  "On  the  Validity 
of  Continuum  Theory  for  Satellite  and  Hypersonic  Flight 
Problems  at  High  Altitudes,"  Jet  Propulsion.  Vol.  28 
(February  1958),  pp.  86-90. 

3.  The  estimated  properties  of  the  upper  atmosphere  are  undergoing 
considerable  revision  as  a  result  of  satellite  measurements.  For 
a  discussion  of  some  of  these  measurements  and  references,  see 

La  Gow,  H.  E.  and  Horowitz,  R.  "Comparison  of  High 
Altitude  Rocket  and  Satellite  Density  Measurements," 

Phys.  of  Fluids.  Vol.  1  ( November-December  1958), 

pp.  AiIFATT. 

4.  A  few  measurements  of  the  variation  with  speed  of  the  accommodation 
coefficients  of  several  gases  on  engineering  surfaces  are  given  in 

Devienne,  F.  M,  "An  Experimental  Study  of  the  Stag¬ 
nation  Temperature  in  a  Free  Molecular  Flow,"  J. 

Aeronaut.  Sci. .  Vol.  24  (June  1957),  pp.  403-4TJ6. 

Rotating  cylinder  measurements  of  the  slip  coefficient  a  indicate 
that  its  value  may  depend  on  the  pressure  "history"  of  the  surface. 
See 


Merlic,  E.  An  Experimental  Determination  of  Reflection 
Coefficients  for  Air  on  Aluminum.  Report  HE-150-141. 
University  of  California  Engineering  Project,  August 
1956. 
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5.  The  aerodynamic  stability  characteristics  of  a  wedge  at  angle  of 
attack  in  free  molecular  flow  have  been  computed.  Results  are 
given  in 

King,  H.  H.  Static  Stability  Analysis  for  a  Wedge 
in  Free  Molecule  Flow.  Report  HE-15Q-165.  Univer¬ 
sity  of  California  Engineering  Project,  December 
1958. 

6.  The  most  recent  experiments  and  theory  indicate  that  the  Navier- 
Stokes  and  Burnett  equations  give  about  the  same  results  in  the 
regime  of  slight  rarefaction.  See 

Talbot,  L.  and  Sherman,  F.  S.  Structure  of  Weak 
Shock  Waves  in  a  Monatomic  Gas.  NASA  Memo.  12-14-58W. 

January  1959. 

Greenspan,  M.  "Propagation  of  Sound  in  Five  Mon¬ 
atomic  Gases,"  J.  Acoust.  Soc.  Am..  Vol.  28  (July 
1956),  pp.  644-S48. 

7.  Local  heat  transfer  measurements  on  the  surface  of  a  circular  cy¬ 
linder  perpendicular  to  a  low  density  flow  are  reported  in 

Tewfik,  0.  and  Giedt,  W.  H.  Heat  Transfer,  Recovery 
Factor  and  Pressure  Distributions  Around  a  Cylinder 
Normal  to  a  Supersonic  Rarefied  Gas  Stream.  Report 
HE-15Q-162.  University  of  California  Engineering 
Project,  January  1959. 

8.  Considerable  data  are  available  for  the  induced  pressures  on  flat 
plates.  See,  for  example, 

Aroesty,  J.  Induced  Pressures  on  Flat  Plates  at 
M  =  4  in  Low  Density  Flow.  Report  HE-15Q-156. 

University  of  California  Engineering  Project,  July 
1958. 

9.  Pressure  distributions  on  spherically-blunted  cones  in  low  density 
flow  are  reported  in 

Talbot,  L.,  Schaaf,  S.  A.,  and  Hurlbut,  F.  C.  "Pres¬ 
sure  Distributions  on  Blunt-Nosed  Cones  in  Low  Density 
Hypersonic  Flow."  Jet  Propulsion,  Vol.  28  (December 
1958),  pp.  832-833. 

Aerodynamic  characteristics  of  a  cone-cylinder  configuration  are 
given  in 

Nark,  T.  Lift  and  Drag  on  Cone  Cylinders.  Report 
HE-150-154.  University  of  California  Engineering 
Project,  December  1957. 
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Lift,  drag,  and  pitching  moment  characteristics  of  a  pointed  cone 
at  angle  of  attack  in  a  low  density  flow  are  reported  in 

Langelo,  V.  and  Lengyel,  A.  Aerodynamic  Coeffi¬ 
cients  of  a  9°  Half-Angle  Pointed  Cone  in  Rarefied 
Gas  Flow,  Report  No.  538-270.  General  Electric 
Company,  May  1958. 

10.  Additional  data  on  impact  probe  corrections  can  be  found  in 

Matthews,  M.  L.  An  Experimental  Investigation  of 
Viscous  Effects  on  Static  and  Impact  Pressure  Probes 
in  Hypersonic  Flow.  GALCIT  Hypersonic  Research  Project 
Memorandum  No.  44.  Guggenheim  Aeronautical  Laboratory, 
California  Institute  of  Technology,  June  1958. 
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• 

1601.3  Table  I 

Thermal  accommodation  coefficient,  a , 

for  air  (Ref.  7) 

Surface 

5 

Flat  lacquer  on  bronze 

0. 88-0.89 

Polished  bronze 

0.91-0.94 

Machined  bronze 

0.89-0.93 

Etched  bronze 

0.93-0.95 

• 

Polished  cast  iron 

0.87-0.93 

Machined  cast  iron 

0.87-0.88 

Etched  cast  iron 

0.89-0.96 

Polished  aluminum 

0.87-0.95 

Machined  aluminum 

0.95-0.97 

Etched  aluminum 

0.89-0.97 

1601.3  Table  II 

Values  of  reflection  coefficient. 

a,  (Ref.  8) 

Gas  and  Surface  Condition 

a 

• 

Air  or  CC^  on  machined  brass  or  old  shellac  1.00 

Air  on  oil 

0.895 

CC>2  on  oil 

0.92 

H2  on  oil 

0.93 

Air  on  glass 

0.89 

He  on  oil 

0.87 

Air  on  fresh  shellac 

0.79 

• 

Altitude 
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1601.2  Figure  1.  The  regimes  of  gas  dynamics. 
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1602.31  Figure  1.  Modified  recovery  factor  for  the  front  side 
of  an  insulated  flat  plate  in  free  molecular  flow. 
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1602.31  Figure  2.  Modified  recovery 
of  an  insulated  flat  plate  in  f 
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1602,31  Figure  3.  Modified  recovery  factor  for  a  flat  plate  in  free  molecular 
flow  with  the  front  and  back  surfaces  in  thermal  contact. 
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1602.31  Fig.  4 
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1602.31  Figure  4.  Modified  Stanton  number  for  the  front  side 
of  an  insulated  flat  plate  in  free  molecular  flow. 
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1602.31  Figure  5.  Modified  Stanton  number  for  the  bac 
of  an  insulated  flat  plate  in  free  molecular  flov 
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1602.31  Figure  6.  Modified  Stanton  number  for  a  flat  plate  in  free  molecular 
flow  with  the  front  and  back  surfaces  in  thermal  contact. 
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1602.32  Figure  1.  Modified  recovery  factor  for  a  sphere  and  a  cylinder 
transverse  to  the  flow  direction  in  free  molecular  flow. 
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1602.32  Fig.  2 


1602.32  Figure  2.  Modified  Stanton  number  for  a  sphere  and  a  cylinder 
transverse  to  the  flow  direction  in  free  molecular  flow. 
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1602.32  Figure  3.  Functions  for  the  determination  of  the  recovery  factor  and  the 
Stanton  number  for  the  cylinder  as  a  function  of  the  angle  of  attack,  a  . 
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1602.4  Figure  1.  Lift  coefficients  for  a  flat  plate 
in  free  molecular  flow:  <y  -  <r'  =0. 
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1602.4  Figure  2.  Lift  coefficients  for  a  flat  plate 
in  free  molecular  flow:  cr  =  cr’  =1  and  T  =  T^. 


1602.4  Figure  3.  Lift  coefficients  for  a  flat  plate  in  free  molecular 
flow:  cr  =  a’  =  1  and  Tw  =  Tw{eq>). 
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1602.4  Figure  4.  Drag  coefficients  for  a  flat  plate  in  free 
molecular  flow:  c 7  =<r’  =  0. 
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1602.4  Figure  5.  Drag  coefficients  for  a  flat  plate  in  free 
molecular  flow:  a  =  a*  =1  and  Tw  =  . 


1602.4  Figure  7.  Drag  coefficients  for  a  sphere  and  a  cylinder 
transverse  to  the  flow  direction  in  free  molecular  flow. 
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1602.4  Figure  8.  Lift  coefficient  for  a  cone  in 
free  molecular  flow:  <r  =  <r'  =0  (Ref.  15). 
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1602.4  Figure  10.  Drag  coefficient  for  a  cone  in  free 
molecular  flow:  a  =  c1  =0  (Ref.  15). 


1602.4  Figure  11.  Drag  coefficient  for  a  cone  in  free 
molecular  flow:  cr  =  <r '  =1  and  Tw  =  T,^  (Ref.  15). 
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1603.41  Figure  2.  Average  recovery  factors  for  spheres  in  supersonic  flow. 
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1603.42  Figure  2.  Variation  of  equilibrium  temperature  for  right  cir 
cylinders  transverse  to  the  flow.  Taken  from  experimental  data 
(Ref.  29) .  M  =  1.3  to  4.5. 
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1603.43  Figure  2.  Thermal  recovery  factors  for  cones  in 

supersonic  flow. 
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1603.51  Figure  1.  Drag  coefficient  for  a  sphere 
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1603.53  Fig.  1 


1603.53  Figure  1.  Surface  pressures  on  a  cone  at  M~4.  (Note:  The  characteris¬ 
tic  length,  x,  is  measured  from  the  vertex  along  the  cone  surface,  and  and 

(Re  ) 9  are  evaluated  using  ideal  flow  properties  behind  a  conical  shock  wave. 


1603.54  Figure  1.  Base  pressure  coefficient  for  a  cone-cylinder  configuratio 
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1603.55  Figure  1.  Increment  of  measured  pressure  coefficient 
over  ideal  value  in  subsonic  flow  ( 0. 10 < M  < 0. 67) . 
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1603.55  Figure  2.  Correction  factors  for  impact  pressure  measure¬ 
ments  for  source-shaped  tubes  in  supersonic  flow. 
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1603.55  Figure  3.  Correction  factors  for  impact  pressure 
measurements  for  open-ended  tubes  in  supersonic  flow. 
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